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Abstract. Let (V, 5) be a finite metric space, where V is a set of n points and 5 is 
a distance function defined for these points. Assume that (V, 5) has a constant doubling 
dimension d and assume that each point p £ V has a disk of radius r(p) around it. The disk 
graph that corresponds to V and r(-) is a directed graph I(V,E,r), whose vertices are the 
points of V and whose edge set includes a directed edge from p to q if S(p, q) < r(p). In [8] 
we presented an algorithm for constructing a (1 + e)-spanner of size 0(n/e d log M), where 
M is the maximal radius r(p). The current paper presents two results. The first shows that 
the spanner of [8] is essentially optimal, i.e., for metrics of constant doubling dimension it 
is not possible to guarantee a spanner whose size is independent of M. The second result 
shows that by slightly relaxing the requirements and allowing a small perturbation of the 
radius assignment, considerably better spanners can be constructed. In particular, we show 
that if it is allowed to use edges of the disk graph I(V, E, ri+ E ), where n+ e (p) = (l + e)-r(p) 
for every p £ V, then it is possible to get a (1 + e)-spanner of size 0(n/e ) for I(V, E, r). 
Our algorithm is simple and can be implemented efficiently. 
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Introduction 

This paper concerns efficient constructions of spanners for disk graphs, an important 
family of directed graphs. A spanner is essentially a skeleton of the graph, namely, a sparse 
spanning subgraph that faithfully represents distances. Formally, a subgraph H of a graph 
G is a t-spanner of G if 5h(u,v) < t ■ 5c(u,v) for every two nodes u and v, where 5g'(u,v) 
denotes the distance between u and v in G' . We refer to t as the stretch factor of the spanner. 
Graph spanners have received considerable attention over the last two decades, and were 
used implicitly or explicitly as key ingredients of various distributed applications. It is 
known how to efficiently construct a (2k — l)-spanner of size 

0( n i+V*) for every weighted 

undirected graph, and this size-stretch tradeoff is conjectured to be tight. Baswana and 
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Sen [?] presented a linear time randomized algorithm for computing such a spanner. In 
directed graphs, however, the situation is different. No such general size-stretch tradeoff 
can exist, as indicated by considering the example of a directed bipartite graph G in which 
all the edges are directed from one side to the other; clearly, the only spanner of G is G 
itself, as any spanner for G must contain every edge. 

The main difference between undirected and directed graphs is that in undirected graphs 
the distances are symmetric, that is, a path of a certain length from utov can be used also 
from v to u. In directed graphs, however, the existence of a path from u to v does not imply 
anything on the distance in the opposite direction from v to u. Hence, in order to obtain 
a spanner for a directed graph one must impose some restriction either on the graph or on 
its distances. In order to bypass the problem of asymmetric distances of directed graphs, 
Cowen and Wagner [5] introduced the notion of roundtrip distances in which the distance 
between u and v is composed of the shortest path from u to v plus the shortest path from 
v to u. It is easy to see that under this definition distances are symmetric also in directed 
graphs. It is shown by Cowen and Wagner [5] and later by Roditty, Thorup and Zwick [6] 
that methods of path approximations from undirected graphs can work using more ideas 
also in directed graphs when roundtrip distances are considered. Bollobas, Coppersmith 
and Elkin [?] introduced the notion of distance preservers and showed that they exist also 
in directed graphs. 

In [8] we presented a spanner construction for directed graphs without symmetric dis- 
tances. The restriction that we imposed on the graph was that it must be a disk graph. 
More formally, let (V, 5) be a finite metric space of constant doubling dimension d, where 
V is a set of n points and 5 is a distance function defined for these points. A metric is said 
to be of constant doubling dimension if a ball with radius r can be covered by at most a 
constant number of balls of radius r/2. Every point p £ V is assigned with a radius r{p). 
The disk graph that corresponds to V and r(-) is a directed graph I(V, E, r), whose vertices 
are the points of V and whose edge set includes a directed edge from p to q if q is inside 
the disk of p, that is, 8{p,q) < r(p). In [8] we presented an algorithm for constructing a 
(1 + e)-spanner with size 0(n/e d log M), where M is the maximal radius. In the case that 
we remove the radius restriction the resulted graph is the complete undirected graph where 
the weight of every edge is the distance between its endpoint. In such a case it is possible 
to create (1 + e)-spanners of size 0{n/e d ), see 0], [2] and [9] for more details. Moreover, 
when the radii are all the same and the graph is the unit disk graph then it is also possible 
to create (1 + e)-spanners of size 0{n/e d ), see [3], [8]. 

As a result of that, a natural question is whether a spanner size of 0{n/e d \ogM) m 
the case of directed disk graph is indeed the best possible or maybe it is possible to get 
a spanner of size 0(n/e d ) as in the cases of the complete graph and the unit disk graph. 
For the case of the Euclidean metric space, the answer turns out to be positive; a simple 
modification of the Yao graph construction fJT] to fit the directed case yields a directed 
spanner of size 0{n/e d ). However, the question remains for more general metric spaces, 
and in particular for the important family of metric spaces of bounded doubling dimension. 

In this paper we provide an answer for this question. We show that our construction 
from [8] is essentially optimal by providing a metric space with a constant doubling dimen- 
sion and a radius assignment whose corresponding disk graph has f2(n 2 ) edges and none of 
its edges can be removed. (This does not contradict our spanner construction from [8] as 
the maximal radius in that case is 0(2 n ) and hence logM = n.) 
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This (essentially negative) optimality result motivates our main interest in the current 
paper, which focuses on attempts to slightly relax the assumptions of the model, in order 
to obtain sparser spanner constructions. Indeed, it turns out that such sparser spanner 
constructions are feasible under a suitably relaxed model. Specifically, we demonstrate the 
fact that if a small perturbation of the radius assignment is allowed, then a (1 + e)-spanner 
of size 0(n/e d ) is attainable. More formally, we show that if we are allowed to use edges 
of the disk graph I(V, E,r± +e ), where r\ +e (p) = (1 + e) • r{p) for every p £ V, then it is 
possible to get a (1 + e)-spanner of size 0(n/e d ) for the original disk graph I(V, E, r). This 
approach is similar in its nature to the notation of emulators introduced by Dor, Halperin 
and Zwick [1] . An emulator of a graph may use any edge that does not exist in the graph in 
order to approximate its distances. It was used in the context of spanners with an additive 
stretch. 

The main application of disk graph spanners is for topology control in the wireless 
ad hoc network model. In this model the power required for transmitting from p to q is 
commonly taken to be 5(p, q) a , where 5(p, q) denotes the distance between p and q and a is 
a constant typically assumed to be between 2 and 4. Most of the ad hoc network literature 
makes the assumption that the transmission range of all nodes is identical, and consequently 
represents the network by a unit disk graph (UDG), namely, a graph in which two nodes 
p, q are adjacent if their distance satisfies 5(p, q) < 1. A unit disk graph can have as many 
as 0(n 2 ) edges. 

There is an extensive body of literature on spanners of unit disk graphs. Gao et al. [3], 
Wang and Yang-Li [10] and Yang-Li et al. [7] considered the restricted Delaunay graph, 
whose worst-case stretch is constant (larger than 1 + e). In [8] we showed that any (1 + e)- 
geometric spanner can be turned into a (1 + e)-UDG spanner. 

Disk graphs are a natural generalization of unit disk graphs, that provide an intermedi- 
ate model between the complete graph and the unit disk graph. Our size efficient spanner 
construction for disk graphs whose radii are allowed to be slightly larger falls exactly into 
the model of networks in which the stations can change their transmission power. In partic- 
ular our constriction implies that if any station increases its transmission power by a small 
fraction then a considerably improved topology can be built for the network. 

Our result has both practical and theoretical implications. From a practical point of 
view it shows that, in certain scenarios, extending the transmission radii even by a small 
factor can significantly improve the overall quality of the network topology. The result 
is also very intriguing from a theoretical standpoint, as to the best of our knowledge, our 
relaxed spanner is the first example of a spanner construction for directed graphs that enjoys 
the same properties as the best constructions for undirected graphs. (As mentioned above, 
it is easy to see that for general directed graphs, it is not possible to have an algorithm 
that given any directed graph produces a sparse spanner for it.) In that sense, our result 
can be viewed as a significant step towards gaining a better understanding for some of the 
fundamental differences between directed and undirected graphs. Our result also opens 
several new research directions in the relaxed model of disk graphs. The most obvious 
research questions that arise are whether it is possible to obtain other objects that are 
known to exist in undirected graphs, such as compact routing schemes and distance oracles, 
for disk graphs as well. 

The rest of this paper is organized as follows. In the next section we present a metric 
space of constant doubling dimension in which no edge can be removed from its corre- 
sponding disk graph. Section [2] first describes a simple variant of our construction from [8], 
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Figure 1: (a) First step in constructing the non-sparsifiable disk graph G. (b) The non- 
sparsifiable disk graph G. 



and then uses it together with new ideas in order to obtain our new relaxed construction. 
Finally, in Section [3] we present some concluding remarks and open problems. 



1. Optimality of the spanner construction 

In this section we build a disk graph G with 2n vertices and Q(n 2 ) edges that is non- 
sparsifiable, namely, whose only spanner is G itself. In this graph M = Q(2 n ) hence our 
spanner construction from [8] has a size of Q(n 2 ) and is essentially optimal. 

Given a set of points, we present a distance function such that for a given assignment 
of radii for the points any spanner of the resulting disk graph must have fi(n 2 ) edges. We 
then prove that the underlying metric space has a constant doubling dimension. 

We partition the points into two types, Y = {y±, . . . , y n } and X = {xi, . . . , x n }. We 
now define the distance function #(■,•) and the radii assignment r(-). The main idea is to 
create a bipartite graph G(X, Y, E) in which every point of Y is connected by a directed 
edge to all the points of X. 

The distance between any two points Xi and Xj is at least 1 + e for some small < e < 1 
and the radius assignment of every point Xi is exactly 1. Thus, there are no edges between 
the points of X. 

We now define the distances between the points of Y and the points of X. We start 
with the point y%. Let 6(yi,Xi) = n for every X{ G X and let r(y%) = n. Place the points of 
X on the boundary of a ball of radius n centered at y\ such that the distance between any 
two consecutive points X{ and Xi+i is exactly 1 + e. This is depicted in Figure Ufa). 

Turning to the point g/2> let 5(y2,Xi) = 2n for every Xi £ X, 5(y2,yi) = 2n + e, and 
r {lj2) = 2n. Hence there is an edge from yi to all the points of X, but no edge connects y2 
and y%. 

We now turn to define the general case. Consider j/j € Y. Let r(yi) = 2 i_1 n and 
S(yi,Xj) = 2*~ 1 n for every Xj G X. Let S(yi,yi-i) = 2 t ~ 1 n + e, and in general, for every 
< j < i we have 

S(yi,yj) = ^2$(yk+i,yk) , (l-i) 

k=j 
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implying that 

5( yi , yj ) < 2*n. (1.2) 
It is easy to verify that yi has outgoing edges to the points of X (and to them only) and it 
does not have any incoming edges. See Figure H{b). 

The resulting disk graph G has 2n vertices and Q(n 2 ) edges. Clearly, removing any 
edge from G will increase the distance between its head and its tail to infinity, and thus the 
only spanner of G is G itself. 

It is left to show that the metric space defined above for G has a constant doubling 
dimension. Given a metric space (V, 5), its doubling dimension is defined to be the minimal 
value d such that every ball B of radius r in the metric space can be covered by 2 d balls of 
radius r/2. In the next Theorem we prove that for the metric space described above, d is 
constant. 

Theorem 1.1. The metric space (XUY,5) defined for G has a constant doubling dimension. 

Proof. Let B be a ball with an arbitrary radius r. We show that it is possible to cover all 
the points of X U Y within B using a constant number of balls whose radius is r/2. The 
proof is divided into two cases. 

Case a: There is some yj G Y within the ball B. (If there is more than one such 
point, then let yj be the point whose index is maximal.) Let B' be a ball of radius R = 2r 
centered at yj. Clearly B C B', so B' contains all the points of B. In what follows we 
show that all the points of X U Y within B' can be covered by a constant number of balls 
of radius r/2. Let y\ be the point within B' whose index is maximal. We have to consider 
two possible scenarios. The first is that yj = y^. This implies that 2/7-4. i ^ B', hence 
R < 5(yj + i,yj) = 2 J n + e. We now show that it is possible to cover B' by a constant 
number of balls of radius R/4. If R < 2 J ~ 1 n, then only yj is within B' and it is covered by 
a ball of radius R/4 centered at itself. If 2- ?_1 n < R < 2 J ~ 1 n + e, then B' contains all the 
points of X and yj. From packing arguments it follows that it is possible to cover all the 
points of X by a constant number of balls of radius n/4, hence also by a constant number 
of balls of radius R > n. The point yj itself is covered by a ball centered at it. Finally, if 
2- ?_1 n + e < R < 2 J n + e, then R/4 is at least 2 J_3 n + e/4. A ball centered at 2/7-3 of radius 
R/4 covers every y^ within B', where 1 < k < j — 3, as <5 (2/7—3,2/*;) < 2 J_3 n. Hence, we cover 
Y n B' by balls of radius R/4 whose centers are yj, yj-i, yj-2 and 2/j-3- We cover X n B' 
as before. This completes the first scenario, where yi = yj. Assume now that yi ^ yj. This 
implies that 5(yt,yj) < R and that R < 5(yi+i,yj), where the first inequality follows from 
the fact that yi € B' and the second inequality follows from the fact that j/i is the point 
with maximal index inside B', hence, y%+\ ^ B' . As 6(yi,yi-i) < 6(yi,yj), we get that 
2*-i n + e < R. Also, by ljT2]> . S(y i+1 , yj ) < 2 i+1 n. We conclude that 2 i ~ 1 n < R < 2 i+1 n 
and that R/4 > 2 l ~ 3 n. A ball centered at yi-3 of radius R/4 covers every y^. within B', 
where k < i — 3, as <5(yi-3,?/fe) < 2 l_3 n. Hence, we can cover B' D Y by balls of radius R/4 
whose centers are j/j, j/j-i, y%-2 and yi-3. We cover X n B' as before. This completes the 
first case. 

Case b: The ball B does not contain any point from Y. The points of X are spread 
as appears in Figure [Ha), thus by standard packing arguments, any ball that contains only 
points from X is covered by a constant number of balls of half the radius. ■ 
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2. Improved spanner in the relaxed disk graph model 

The (negative) optimality result from the previous section motivates us to look for a 
slightly relaxed definition of disk graphs in which it will still be possible to create a spanner 
of size 0{n/e d ). 

Let (V, 5) be a metric space of constant doubling dimension d with a radius assignment 
r(-) for its points and let / = (V,E,r) be its corresponding disk graph. Assume that we 
multiply the radius assignment of every point by a factor of 1 + e, for some e > 0, and 
let I' = (V,E',ri +e ) be the corresponding disk graph. It is easy to see that E C E'. In 
this section we show that it is possible to create a (1 + e)-spanner of size 0(n/e d ) if we are 
allowed to use edges of I'. As a first step we present a simple variant of our (1 + e)-spanner 
construction of size 0(n/e d log M) from [8j. This variation is needed in order to obtain the 
efficient construction in the relaxed model which is presented right afterwards. 

2.1. Spanners for general disk graphs 

Let (V, 5) be a metric space of constant doubling dimension and assume that any point 
p G V is the center of a ball of radius r(p), where r{p) is taken from the range [1,M]. In 
this section we describe a simple variant of our construction from [8], which computes a 
(1 + e)-spanner with 0(n/e d log M) edges for a given disk graph. We then use this variant, 
together with new ideas, in order to obtain (in the next section) our main result, namely, a 
spanner with only 0(n/e d ) edges. 

The spanner construction algorithm receives as input a directed graph I(V, E, r) and 
an arbitrarily small (constant) approximation factor e > 0, and constructs a set of span- 
ner edges Egp R , returning the spanner subgraph H mR (V, Egp R ). The construction of the 
spanner is based on a hierarchical partition of the points of V that takes into account the 
different radius of each point. The construction operates as follows. Let a and f3 be two 
small constants depending on e, to be fixed later on. Assume that the ball radii are scaled 
so that the smallest edge in the disk graph is of weight 1. Let i be an integer from the 
range [0, |d°gi+ a M\\ and let M, = Mj (1 + cif . The edges of I(V, E, r) are partitioned into 
classes by length, letting E(M i+ i,Mi) = {(x,y) \ M i+ i < 5(x,y) < Mi}. Let £(x,y) be the 
level of the edge (x, y), that is, £(x, y) = i such that (x, y) G E(M,i + \, Mi). Let p be a point 
whose ball is of radius r{p) G [Mj + i,Mj]. It follows that level i is the first level in which p 
can have outgoing edges. We denote this level by £(p). 

For every i G [0, [log 1+Q ,Mj], starting from i = 0, the edges of the class E(Mi + i, Mi) 
are considered by the algorithm in a non-decreasing order. (Assume that in each class the 
edges are sorted by their weight.) In each stage of the construction we maintain a set of 
pivots Pi. Let x £ V and let NN(x,i-i) be the nearest neighbor of x among the points of 
Pi. For a pivot p G Pi, define Ti{p) = {x \ x G V, NN(x, Pi) = p, r(x) > 5(x,p)}, namely, all 
the points that have a directed edge to p and p is their nearest neighbor from Pi. We refer 
to Ti(p) as the close neighborhood of p. 

The algorithm is given in Figure [2j Let (x,y) be an edge considered by the algorithm 
in the ith iteration. The algorithm first checks whether x or y or both should be added 
to the pivots set Pj. The main change with respect to [8] is that if y is assigned with 
a large enough radius it might become a pivot when the edge (x, y) is examined. When 
considering the edge (x, y), the algorithm acts according to the following rule: If the distance 
from x to its nearest neighbor in Pi is greater than f3Mi + \ then x is added to Pi. If the 
distance from y to its nearest neighbor in Pi is greater than /3Mj + i and the radius of y 
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Algorithm disk-spanner (I(V,E,R),e) 

£ S D P R <- 4> 

for i <- to [log 1+Q Mj 

for each (x,y) G P(Mj + i,Mj) do 

if S(NN(x,Pi),x) > (3M i+1 then P <- P U {x} 

if <5(NN(y, P;),y) > /3M m A r(y) > M i+1 then P <- P U {y} 

if r(y) > M m 

if ftsW) € -Eg ) p R s.t. x' G r,(NN(x,P,))A y' G r,(NN(y,P)) 

then Egp R <- P™ R U {(x, y)} 
if r(y) < M i+ i 
if J(x',y) G PgP s.t. x 1 G r i (NN(x,P l )) 

thenP s D P R ^P s D P R U{(x,y)} 

Pi+l ^— Pi 

return P DIR (V, 



Figure 2: A high level implementation of the spanner construction algorithm for general 
disk graphs 



is at least Mj+i then y is added to Pj. To decide whether the edge (x,y) is added to the 
spanner, the following two cases are considered. The first case is when r(y) > Mj + i. In 
this case, if there is no edge from the close neighborhood of x to the close neighborhood 
of y then (x,y) is added to the spanner. The second case is when r(y) < Mj+i. In this 
case, if there is no edge from the close neighborhood of x to y then (x, y) is added to the 
spanner. When i reaches [logi +Q M\ , the algorithm handles all the edges that belong to 
E(M\\ os a mj+1) -^Liogi+c* M \ )■ This includes also edges whose weight is 1, the minimal 
possible weight. The algorithm returns the directed graph H Bm (V, P°p R ). 

In what follows we prove that for suitably chosen a and f3, H DIR (V, P<Pp R ) is a (1 + e)- 
spanner with 0{n/e d log M) edges of the directed graph I(V,E,r). 

Lemma 2.1 (Stretch). Let e > 0, set a = /3 < e/6 and let H = H DIR (V,E^ R ) be the 
graph returned by Algorithm disk- spanner (P^, E, r), e). // (x,y) G E then <5f/(x,y) < 
(l + e)(5(x,y). 

Proof. Recall that the radii are scaled so that the shortest edge is of weight 1. We prove 
that every directed edge of an arbitrary node x G V is approximated with 1 + e stretch. 
Let i G [0, [log 1+Q MJ]. The proof is by induction on i. For a given node x, the base of the 
induction is the maximal value of i in which x has an edge in P(Mj + i, Mi). Let j be this 
value for x, that is, the set E(Mj + ±, Mj) contains the shortest edge that touches x. Every 
other node is at distance at least Mj+i away from x, hence x is a pivot at this stage and 
every edge that touches x from the set E(Mj + i,Mj) is added to E®p R . 

Let (x,y) G P(Mj + i,Mj) for some i < j and let p = NN(x,Pj). Assume that r(y) > 
Mi + \ and let q = NN(y,Pj). It follows from definition that S(x,p) < j3Mi + \ and 5(y,q) < 
PM i+1 . 

If the edge (x, y) is not in the spanner, then there must be an edge (x, y) G Pgp R , where 
x G Tj (p) and y G Tj (q) . The crucial observation is that the radius of x and y is at least 
M i+1 . By the choice of /3, it follows that 2(3M i+ i < M i+ i and (x,x),(y,y) G E. Thus, 
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there is a (directed) path from x to y of the form (x, x, y, y) whose length is 4/3Mj+i + Mj. 
However, only its middle edge, (x, y), is in Efp R . The length of this edge is bounded by the 
length of the edge (x, y) since the algorithm picked the minimal edge that connects between 
the neighborhoods. This implies that the length of (x,y) is at most Mj. 

By the inductive hypothesis, the edges (x, x) and (y, y) whose weight is at most 2/3Mj + i 
are approximated with 1 + e stretch. Thus, there is a path in the spanner from x to y whose 
length is at most (1 + e)5(x, x) + Mi + (1 + e)5(y, y), and this can be bounded by 

(1 + e)4/3M i+1 + Mi = ((l + e)4/3 + (l + a))M +1 . 

As the edge (x,y) G E(Mi + i, Mi) it follows that 5(x,y) > Mj+i. It remains to prove that 
1 + 4e/3 + 4/3 + a < 1 + e, which follows directly from the choice of a and /3. 

If r(y) < Mj + i then there must be an edge (x,y) G E$p R , where x G Ti{p). Following 
similar arguments to those used above it can be shown that there is a path in the spanner 
from x to y of length at most (1 + e)2/3Mj + i + Mi and hence bounded by (1 + e)Mi + \. m 

The size of the spanner. We now prove that the size of the spanner H mR (V, E$p R ) is 
0(n/e d log M). As a first step, we state the following well-known lemma, cf. [2]. 

Lemma 2.2. [Packing Lemma] If all points in a set U G M. d are at least r apart from each 
other, then there are at most (2R/r + l) d points in U within any ball X of radius R. 

The next lemma establishes a bound on the number of incoming spanner edges that a 
point may be assigned on stage i G [0, Llog 1+Q , MJ] of the algorithm. 

Lemma 2.3. Let i G [0, Ll°gi+a MJ] and let y G V. The total number of incoming edges of 
y that were added to the spanner on stage i is 0(e~ d ). 

Proof. Let (x,y) be a spanner edge and let NN(x,Pi) = p. We associate (x,y) to p. From 
the spanner construction algorithm it follows that this is the only incoming edge of y whose 
source is in Ti(p). Thus, this is the only incoming edge of y which is associated to p. 
Now consider all the incoming edges of y on stage i. The source of each of these edges is 
associated to a unique pivot within distance of at most Mj + 2/3Mj + i away from y and any 
two pivots are /3Mj + i apart from each other. Using Lemma 12.21 we get that the number of 
edges entering y is ( M '+g^' +1 + l) d = ((1 + a)/P + 3) d = 0{e~ d ). m 

It follows from the above lemma that the total number of edges that were added to 
El?p R in the main loop is 0(n/e d log M). The total cost of the construction algorithm is 
0{m log n). For more details on the construction time see [8]. 

2.2. Spanner for relaxed disk graphs 

Let (V, 5) be a metric space of constant doubling dimension d with a radius assignment 
r(-) for its points and let I = (V,E,r) be its corresponding disk graph. Assume that we 
multiply the radius assignment of every point by a factor of 1 + e, for some e > 0, and 
let /' = (V, E' ,ri +t ) be the corresponding disk graph. In this section we show that it is 
possible to create a (1 + e)-spanner of I of size 0(n/e d ) if we are allowed to use edges of I'. 

Our construction consists of two stages: a building stage and a pruning stage. The 
building stage creates two spanners, H and H' , using the algorithm of Section [2.11 where 
H is the spanner of I and H' is the spanner of In the pruning stage we prune the union 
of these two spanners. Throughout the pruning stage we use the radius assignment of each 
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point before the increase. Let q € V and let £(q) be the first level in which q can have 
outgoing edges, that is, r{q) £ [M^ g ) +1 , M^i q \] (recall that as the levels get larger the edges 
get shorter). In the pruning stage we only prune incoming edges of q whose level is below 
£{q). In other words, we do not touch the incoming edges of q that are shorter than the 
radius of q. The pruning is done as follow. Let 7 = log 1+a 1/(3 + 1. We keep in the spanner 
the incoming edges of q that come from the first 47 different levels below £{q). 

Let H be the resulting spanner and let E be the remaining set of edges after the pruning 
step. In the remainder of this section we show that the size of H is 0{n/e d ) and its stretch 
with respect to the distances in /(V, E, r) is 1 + e. We start by showing that the size of H is 
0(n/e d ). Notice that the first part of the proof below is possible only due to the change we 
have done in the previous section to our spanner construction from [8]. Roughly speaking, 
given an edge (p, q) € E that is shorter than r(g) we use pivot selection also on g's side 
(and not only on p's) to sparisify the graph. This allows us to deal separately with edges 
of q of length larger than r(q) and those of length smaller than r(q). 

Lemma 2.4. \E\ = 0(n/e d ). 

Proof. Let (p, q) be a spanner edge that survived the pruning step. There are two possible 
cases to consider. 

The first case is that £(p,q) > £(q). Let i = £(p,q) and let x = NN(p, Pi) and y = 
NN(g, Pi). By packing considerations similar to Lemma [2.3l it follows that the total number 
of edges at level i that connects between two pivots as the edge (p, q) that are associated 
with x (and with y) is 0(l/e d ). The distance between x and y is at most 2/3Mj + i + Mj, 
therefore at level i — 2j either x or y are no longer pivots. 

Let x £ Pj and x Pj-i, that is, Pj is the first pivot set that contains x. Then we 
charge x with every (incoming and outgoing) edge of this type from levels [j, j + 27] that 
is incident to x. Now given such an edge (p, q) whose level is i, either x or y are not pivots 
in level i — 27, which means that either x or y has been charged for this edge, since one of 
them first becomes a pivot between levels i — 2j and i. 

The second case is that £(p, q) < £{q). In this case, it must be that level £(p, q) is among 
the 47 first different levels below £(q) from which an incoming edge is allowed to enter q. 
Subsequently, we associate the edge (p, q) with q, as the total number of such edges that q 
can have is 0{^/e d ). m 

We now turn to prove that the stretch of the spanner H with respect to the disk graph 
/ is 1 + e. 

Lemma 2.5. Let (p,q) be an edge of the spanner H that was pruned. We show that there 
is a path in H whose length is at most (1 + e)5(p, q). 

Proof. The proof is by induction on the lengths of the pruned edges. For the induction base 
let (p, q) be the shortest edge that was pruned. For every x G V, let s(x) be the head of 
an edge whose level is the 7-th level below £(x) from which x has an incoming edge. Let 
qi, . . . qi, . . . be a sequence of points, where qi = q and qi = s(gj_i). As qi + i = s(qi), it 
follows that £(qi + i,qi) < £{qi) — 7. Combining this with the fact that £(qi) < £(qi,qi-i) we 
get that £(q i+ i,qi) < £(q i ,q i - 1 ) - 7. Therefore, 5(q i ,q i - 1 ) < 
The analysis distinguishes between two cases. 

Case a: There is a point qt such that 8(qt,q) > f35(p,q). This situation is depicted 
in Figure El (If there is more than one point that satisfies this requirement, take the one 
whose index is minimal.) 
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Figure 3: The case in which qt exists 
Claim: 5{q u q t -\) > %S(p,q). 

Proof. For the sake of contradiction, assume that 5(qt,qt-i) < ^5{p,q). This implies that 

t 

2% i ,gt-i) < P6(p,q) < S(q t ,q) < J^tffai, , (2.1) 

i=2 

where the last inequality follows from the triangle inequality as the distance between q and 
q t is at most Y,l=2^(.9i-i,Qi)- For every 2 < i < t - 1 we have 8(q i ,q i - 1 ) < P6(q i+ i,qi), 
which implies that S(qi,qi-i) < j3 t ~ l 5{qt,qt~i)- Combined with (|2.1|) . we get 

t-i t-i 
5{q u q t -i) < ^%i,%_i) < 5(gt,gt-i)5^j9*~* • 

i=2 i=2 

If /3 < 1/2 we have Yll=2 ^ ^ ana - ^ ms yi e ld s a contradiction. ■ 

We now focus our attention on the point qt—x- The minimality of qt implies that 
5(q,qt~i) < /36(p,q). By combining it with the triangle inequality we get that 5(p,qt-i) < 
6(p, q) + /3S(p, q). Therefore, in the graph I' there must be an edge from p to qt-i- 

Let i = £(p, qt-i)- There are two possible scenarios for the spanner H'. The first scenario 
is when r'(qt-i) < Mi + \. In this case, there is an edge in H' from some x € Fj(NN(p, i)) to 
qt-i, whose length is at most 6{p,q) + /3S(p,q). 

There are 47 different levels below £(qt-i) from which edges that belong to the spanners 
H and H' are not being pruned and survived to the spanner H. We know that the edge 
(qt,qt-i) is such an edge from the 7-th non-empty level below £(qt~\). We also know that 
5(qt,qt-i) > §5(p, g). Therefore, as the length of the edge (x, qt-i) is at most 5(p,q) + 
P6(p, q) it is within the 47 non-empty levels below £(q t -i) and it is not pruned. We can now 
build a path from p to q by concatenating three segments as follows: A path from p to x, 
the edge (z, qt-i) and a path from qt-i to q. The point x is at most 2/35(p,q) + 2j3 2 5(p,q) 



RELAXED SPANNERS FOR DIRECTED DISK GRAPHS 



619 



away from p and for the right choice of /3 it is less than 8(p,q)/(l + e), hence the weight 
of every edge on the path that approximates the distance between x and p in H U H' is 
less than 5(p,q), the shortest pruned edge, and the entire path survived the punning stage. 
Similarly, the point qt-i is at most f35(p,q) away from q and again for the right choice of 
j3 every edge on the path that approximates the distance between q t ~\ and q survived the 
punning stage. Thus, we get that there is a path whose length is at most 

(1 + e)(3/3S(p, q) + 2(3 2 5(p, q)) + 5(p, q) + 06(p, q) , 

which is less than (1 + e)d(p,q) for f3 < e/11. 

The second scenario is when r'{qt-i) > Mi + \. In this case, there is an edge in H' from 
some x G Tj(NN(p, z)) to some y £ rj(NN(gt_i, i)) whose length is at most S(p, q) + /35(p, q) , 
which is not being pruned. We can build a path from p to q by concatenating three segments 
as follows: A path from p to x, the edge (x, y) and a path from y to q. As before, for the 
right choice of f3 the paths from p to x and from y to q are composed from edges that are 
shorter from S(p,q), the length of the shortest pruned edge, hence, from the minimality 
5(p,q) every edge on these paths survived the punning stage. We get that there is a path 
whose length is at most 

(1 + e)(4/3<f(p, q) + 5/3 2 5(p, q)) + 5(p, q) + P5{p, q) , 

which is less than (1 + e)6(p,q) for (3 < e/19. This completes the proof for case a. 

Case b: There is no point qt such that 5(qt,q) > j35{p,q). In this case, let q t -\ be the 
last point in the sequence of points q\, ., where = s(%_i) and q\ = q. Similarly 
to before, there are two possible scenarios for the spanner H'. Let i = l(p,qt-i). The 
first scenario is when r'(qt-i) < Afj+i. In this case, there is an edge in H' from some 
x G rj(NN(p, i)) to qt-\ whose length is at most 5(p,q) + f36(p,q). This edge could not be 
pruned, since if it was pruned then q t -\ could not have been the last point in the sequence. 
Hence we can construct a path from p to q exactly as we have done in the first scenario 
of case a, described above. The second scenario is when r'(qt-\) > Mi + \. In this case, we 
can construct a path from p to q exactly as we have done in the second scenario of case a, 
described above. 

This completes the proof of the induction base. The proof of the general inductive step 
is almost identical. The only difference is that when a path is constructed from p to q, its 
portions from p to x and from qt-\ to q in the first scenario and from p to x and from y to 
q in the second scenario exist in H by the induction hypothesis and not by the minimality 
oi5{p,q). * 

We end this section by stating its main Theorem. The proof of this Theorem stems 
from Lemma 12.41 and Lemma 12.51 

Theorem 2.6. Let (V, 5) be a metric space of constant doubling dimension with a radius 
assignment r(-) for its points and let I = (V,E,r) be its corresponding disk graph. Let 
I' = (V,E',ri +t ) be the corresponding disk graph in the relaxed model. It is possible to 
create a (1 + e)-spanner of size 0(n/e d ) for I using edges of V . 

3. Concluding remarks and open problems 

This paper presents a spanner construction for disk graphs in a slightly relaxed model 
that is as good as spanners for complete graphs and unit disk graphs. This result opens 
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many other research directions for disk graphs. We list here two questions that we find 
particularly intriguing: Is it possible to design an efficient compact routing scheme for disk 
graphs? And is it possible to build an efficient distance oracle for disk graphs? 
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